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ABSTRACT 

The three flavor version of the Nambu-Jona-Lasinio chiral soliton model for 
baryons is employed to calculate the twist-2 contribution to the polarized 
nucleon structure function gi{x). In particular the role of the strange quark 
degree of freedom as a collective excitation of the chiral soliton is investigated 
in the context of flavor symmetry breaking. The model prediction for gi{x) 
refers to a low momentum scale Qq. The leading order corrections to the 
scale dependence is computed along the QCD evolution program allowing to 
compare with data from SLAG. 
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1. Introduction 

The investigation of the strange quark contribution to the polarized structure function 
is mainly motivated by the empirical results obtained in the context with the proton spin 
puzzle, cf. ref [jl| for a review. In this context one considers nucleon matrix elements of 
axial currents, which are obtained as the zeroth moment of polarized structure functions. 
The puzzle firstly refers to the smallness of the observed nucleon matrix element of the axial 
singlet current. Already early studies in the Skyrme model (the simplest version of a 
chiral soliton model) indicated that chiral soliton models are capable of reproducing that 
result. The data analysis secondly revealed that the strange quark might make a sizable 
contribution to the axial singlet charge of the nucleon; up to a third of that of the down 
quark. This surprising result was obtained using SU(3) symmetric baryon wave-functions; 
the inclusion of symmetry breaking effects reduces this ratio. Certainly it is interesting 
to compute the full dependence of the corresponding structure functions on the Bjorken 
variable. At this point the collective approach to incorporate strange degrees of freedom in 
chiral soliton models P] becomes very attractive^. It not only allows one to account for such 
symmetry breaking effects in the nucleon wave-function but in particular to make explicit 
the strange quark contribution to nucleon structure functions. This is a major advantage 
over other low energy models for the nucleon like e.g. the MIT bag model |^. 

The polarized structure function gi is extracted from the hadronic tensor for electron- 
nucleon scattering, 

= ^Jd"^ ^'"'Hm [MO, 4(0)] \N) , (1) 

where = qiOlfj-Q^liO is the electromagnetic current with Q = diag — |, — being 
the quark charge matrix. The nucleon state, denoted by |A^) in eq (|I|), is characterized by 
its momentum (P^) and spin (S^). The polarized structure function gi parameterizes the 
longitudinal part of the antisymmetric piece, 

Wj^fiQ) = {W^u - W,^)I2% = te^^^^^j^ [g^{x, Q^)S'' + transverse part} . (2) 

Here refers to the momentum transferred to the nucleon by the virtual photon with 
= —q^. The Bjorken variable is defined as x = / P ■ q. We are interested in the leading 
twist contribution to gi. It is extracted from the hadronic tensor by assuming the Bjorken 
limit which corresponds to the kinematic regime 

% = \q\ ~ Mnx with |q| oo and x fixed . (3) 

In this limit eq (0) is straightforwardly inverted yielding 



M 



9i[x) = -rre 



2 ^ (4) 

for the kinematical conditions q = (g'',0,0, \q\), P = (Mat, 0,0,0) and 5* = (0,0,0,1), i.e. 
the nucleon rest frame with the nucleon spin and the transferred momentum being aligned, 
hence the notion longitudinal. 

^ Cf. ref M for a recent review on soliton models in flavor SU(3). 
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In chiral soliton models strange quark effects enter the structure functions in two ways, 
firstly there are direct contributions to the currents which do not exist in the two flavor 
model and secondly the structure of the nucleon state \N) changes when generalizing from 
flavor SU(2) to SU(3). For the present study we will consider the Nambu-Jona-Lasinio 
(NJL) model P] in its bosonized form which is well known to possess a chiral soliton 



solution [|T2] . For comprehensive lists of references we refer to recent reviews [|TT], |T2[ . In this 
model the defining Lagrangian contains only quark fields hence all quantities can formally 
be expressed in terms of these fields. In particular the formal expression for the current is 
that of a free Dirac theory which makes the commutator in eq (|l|) feasible. It is this feature 
which actually allows us to compute structure functions from a chiral soliton. In other 
soliton models, which do not possess such a clear connection to the quark flavor dynamics, 
the computation of structure functions seems infeasible due to the complicated structure 
of the current operator. In the two flavor version of the model various nucleon structure 
functions have been discussed [|T3], |T1|, |T3], 0, 0, [T^ while the present study represents the 



first attempt to consider the effects of strange quarks. 

The extension to flavor SU(3) is particularly interesting because the calculation of 



Qi in the two flavor model already yielded reasonable agreement with experiment. On the 
other hand one expects strange degrees of freedom to have non-negligible impact on nucleon 
properties. 

2. The NJL Chiral Soliton 

In the three dimensional flavor space with up, down and strange quarks the NJL model 
Lagrangian reads 

C = q (z^ - m°) g + 2 Gnjl E ((gy g)' + (gy^75g)') • (5) 

Here q, m° = diag(m°, m^, m°) and Gnjl denote the quark field, the current quark mass 
matrix and a dimensionful coupling constant, respectively. In what follows we will assume 
the isospin limit m° = m^. Of course, we will consider flavor symmetry breaking by allowing 
the strange current quark mass to be different, i.e. m° 7^ m°. 

Using functional integration techniques we obtain the bosonized version of the NJL model 
action ^ 

A[M] = TrA \og{iD) + — ^ / d^x tr |m° (M + M^) - MM^ - {my} , (6) 
D = i^-[m+M^)- 75 {M - M^) . (7) 



The composite scalar (S) and pseudoscalar (V) meson fields are contained in = 5 + iV, 
and appear as quark-antiquark bound states. Apparently A4 represents a 3 x 3 matrix field 
in flavor space which behaves as the sum of scalar {S) and pseudoscalar (V) quark bilinears 
under chiral transformations. For regularization, which is indicated by the cut-off A, we 
will adopt the proper-time scheme |]19|. The free parameters of the model are the current 



quark mass matrix m , the coupling constant Gnjl and the cut-off A. The equation of 
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motion for the scalar field S may be considered as the gap-equation for the order parameter 
(qq) of chiral symmetry breaking. This equation relates the vacuum expectation value, 
(M) = M = diag(M^j,Mrf = Mu^Mg) to the model parameters m°, Gnjl and A. For 
apparent reasons M is called the constituent quark mass matrix. The occurrence of this 
vacuum expectation value reflects the spontaneous breaking of chiral symmetry and causes 
the pseudoscalar fields to emerge as (would-be) Goldstone bosons. At this stage we expand 
A to quadratic order in V around M. Then the model parameters are related to physical 
quantities hke the pion mass, = 135MeV and the pion decay constant, = 93MeV. 
In account of the gap-equation this leaves one undetermined parameter which we choose to 
be the up constituent quark mass Mu 0. Typical values are M„ = 350 ~ 450MeV. The 
kaon mass rriK = 495MeV allows us to fix the strange current quark mass which, subject to 
the gap-equation in the strange sector, also determines the strange constituent quark mass 
Mg. Using the typical values for Mu the model underestimates the kaon decay constant 
/k = 114MeV by about 10-15% [0. 

Turning to the baryon sector of the model we adopt the hedgehog ansatz for the meson 
fields M = ^M^ 

e = Ci/(r-)=expQr-f-e(r)) (8) 

in order to determine the chiral soliton. For static meson configurations as (P) it is straight- 
forward to deduce the classical energy E[Q] functional associated with the action [^| 



+ mlf^ I rfV(l -cose(r)). (9) 

Here refer to the energy eigenvalues of the Dirac Hamiltonian 

ho = a-p + Mupexp(i-f5T-re{r))f + MspS (10) 

which is derived from the operator (|^) D = /3{idt — ho) upon substituting the hedgehog ansatz 
(H). From the appearance of the strange and non-strange projectors S = diagg(0, 0, 1) and 
T = diagf)(l, 1, 0) we observe that the strange quarks are not effected by the hedgehog field. 
In eq (^ the subscript "v" denotes the valence quark level. This state is the distinct level 
bound in the soliton background, i.e. —m < < m. Similar to the energy functional (H) 
other quantities also separate into contributions associated with the explicit occupation of 
the valence level and a (regularized) piece due to the vacuum being polarized by the meson 
fields. The chiral soliton, B(r), is finally obtained by self-consistently extremizing £"[0] 



3. The Nucleon State in Flavor SU(3) 

States with nucleon quantum numbers are generated from the static configuration (^) by 



introducing collective coordinates for the large amplitude fluctuations of the soliton p2|, |[. 
This approach may be viewed as an approximation to the unknown time dependent solution 
to the equations of motion for the meson fields. Subsequently these coordinates are treated 
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quantum-mechanically. The large amplitude fluctuations are associated with the rotations 
in flavor space because the spatial rotations can be absorbed into the former as a consequence 
of the hedgehog ansatz (B). For the two-flavor NJL chiral soliton model this quantization 



approach has been performed in ref [^. In the more involved case of three flavors we have 

m 

ar,t) = A{t)Ur)A^{t) , (11) 

where A{t) is a 3 x 3 matrix in flavor space. Substituting this conflguration and transforming 
to the flavor rotating frame q' = A{t)q reveals that the eigenvalues of the modifled Dirac 
Hamiltonian 



with 



h' = ho + h,ot + hsB (12) 



h,ot = -^yl(t)-^At(t) = l^A„fi" (13) 



hsB = ^ - Ms) rp^Yl DsiX + E ^8"^- + (^88 - 1) As I 

V3 Vi=l a=A ) 



(14) 



enter the functional trace (|^) [^. Besides the angular velocities Of" also the adjoint rep- 
resentation of the collective coordinates Dab = (l/2)tr(AaAAfeA^) have been used to sim- 
plify the additional parts of the Dirac Hamiltonian. For convenience the chiral rotation 
T = {C,Ij + C.h)/'^ + Ini^H ~ ^-H')/2 has been introduced in eq (|I^). Using these deflnitions 
it is straightforward to extract the Lagrangian L{A, f2") for the collective coordinates from 
the action functional (H). As the rotations (0) are assumed to proceed adiabatically the 
expansion of the action functional is terminated at quadratic order in In addition the 
functional trace in (^ is expanded in the difference Mg — which measures the flavor 
symmetry breaking originating from different constituent quark masses^. Essentially this 
represents an expansion in hj-ot + hsB- 

By Legendre transformation to the SU(3) right generators, Rb = —dL{A,Q°')/dQ'' the 
Hamilton operator in the space of the collective coordinates is deduced. It has the form 0] 



H{A,Ra) = 



1 1 



3 8 3 



+^ E Dsi (2/2. + a,Ds.) + A ^ /^^^ (2i?„ + (3,Dsa) + ^7(1- ^ss) 

^'^ 1=1 ■^P a=4 ^ 

+ ^75 (l - Dl,) + E DsiDsi + l^TS E DsaDsa (15) 

together with the constraint Rg = for i? = 1 and Nc = 3. This constraint restricts 

the allowed states to those with half-integer spin. The quantities a^, . . . , 7^5 are functionals 

^The flavor symmetry breaking stemming from different current quark masses resides in tlie local piece 
in (I). 
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of the self-consistent chiral angle. For details of their evaluation and numerical results we 
refer to ref ||2^. Here it is only important to note that for the self-consistent soliton these 



constants of proportionality are dominated by their valence quark contributions. 

The most important feature of the collective Hamiltonian is that it can be diagonalized 
exactly yielding as eigenstates the low-lying and |^ baryons. Due to the presence of 
flavor symmetry breaking these are no longer pure octet (decouplet) states but acquire 
sizable admixture of states with identical quantum numbers in the higher dimensional SU(3) 
representations like 10 and 27. This diagonalization is a generalization of the Yabu-Ando 
approach and is comprehensively described in refs |20, All nucleon matrix elements 



to be computed henceforth will employ these exact eigenstates. 

4- Valence Quark Approximation to gi{x) 

In order to compute gi in the present model we have to bear in mind that we are dealing 
with localized field configurations. The resulting expression for the hadronic tensor in the 
Bjorken limit has been obtained in 0. Its antisymmetric component becomes 

^ e,,u. P sgn (fco) 5 {k') dt S''^^'^-^' 
X j d^xi j d^X2 exp [—i{k + q) ■ {xi — X2)] 

x(Ar|{^(a;i,t)QV7'^(a;2,0) + ^^(a;2,0)QV7'*(a;i,t)}|Ar) , (16) 

where e^pya'^"'^^ is the antisymmetric combination of 7^7p7y. The expression (|16D is es- 
sentially obtained by applying Wick's theorem to the current commutator in eq (|l|). For 
the intermediate quark the free correlator is substituted because in the Bjorken limit this 
quark is far off-shell and hence not sensitive to momenta which are of the scale as those 
attributed to the soliton. Finally a collective coordinate describing the position of the soliton 
in coordinate space is introduced and integrated over. 

We have already noted that the constants of proportionality entering the collective Hamil- 
tonian (|T3p are dominated by their valence quark contribution. This dominance is even more 



pronounced for axial matrix elements like (A^|^7573Aag|iV) JIT], |T2|, For these axial cur- 
rent matrix elements the vacuum contribution is commonly found to be 10% of the total; 
for some flavor combinations even less. This establishes the assumption that the vacuum 
contribution to the polarized structure functions is negligible and manifests itself in the va- 
lence quark approximation for the structure function gi{x). This approximation is defined 
by substituting the valence quark wave-function 

^,{r,t)=e-''^'A{t)^,{r) with <D.(r) = V^.(r) + ^ ^l,,{r) ^^'^-^"^ ^ ^'^'""^ (17) 

in the hadron tensor ([T6|) . The spinors ip^{r) = {r\fi) diagonalize the classical Dirac Hamil- 
tonian ( |T0|) yielding the eigenvalues e^. It is important to stress that the dependence of 



the valence quark wave-function "^^{Vji) on the collective coordinates is included in the 
calculation of gi. Only with this input it will be possible to disentangle the strange quark 
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contribution. The calculation is performed by taking the Fourier transform of \E'v(^; t) which 
allows us to carry out the momentum integrals in eq (0). The technical details of this 
calculation will be presented elsewhere together with results for the unpolarized structure 



functions 24 



The matrix element in eq (0) between nucleon states is to be taken in the space of the 
collective coordinates, A(t) (see eq. ([TT|)) as the object in curly brackets is an operator in 
this space, which can be deduced from eq (|1^. Here it is important to repeat that A(t) 
spans the three dimensional flavor space which brings the strange degrees of freedom into 
the picture. This is the main difference to e.g. bag model calculations where the nucleon 
state is considered to be a product state of three specified quarks. In the present model it is 
rather a collective excitation of quark fields in the background of the self-consistent soliton. 

5. Projection and Evolution 

In the chiral soliton model the baryons states are not momentum eigenstates causing the 
structure functions not to vanish exactly for x > 1. This short-coming is due to the localized 
field configuration and thus the nucleon not being a representation of the Poincare group. 
As a consequence the computed structure functions are not frame independent. It turns out 
that in the infinite momentum frame (IMF) the structure functions indeed vanish for x > 1 



25, 26]. This is a consequence of the Lorentz contraction associated with the boost from 
the rest frame (RF) to the IMF which exactly projects onto < a; < 1. For the polarized 
twist-2 structure function this implies 

^r^(a:) = ^.r(-Mi-x)), (18) 

where Qi^^-' (x) denotes the structure function computed in the nucleon rest frame as discussed 
in the preceding sections. The transformation ([T8|) leaves the integral over gi invariant 
provided that in the rest frame the integration range has been extended to infinity. 

The chiral soliton model is considered to approximate QCD at a low momentum scale Ql 
whence the result gf^^\x) should be interpreted as gi{x, Ql). It should be noted that Ql is 
a new parameter to the model. In order to compare with experimental data it is mandatory 
to evolve this structure function to gi{x, Q^) according to the DGLAP procedure ||2^. Here 

is the momentum scale set by the experiment. In order to apply this procedure we 
first have to separate the flavor singlet (0) and non-singlet {ns) contributions to gi{x,Ql) 
as the former mixes with the gluon contribution g{x,Q'^) upon evolution. This separation 
is straightforwardly accomplished by substituting the appropriate flavor matrix for in 
eq ([TB|). The leading order evolution equations read in the case of three flavors and with 
t = ln(Q^/AQCj-,) as well as oqcd = 47r/9t, 

at 2-71 Jx y \y J 

dgi\x,t) ^ Qqcd dy \ ^ ( ^ ^(o) ^^ , ar> ( 
dt 27r 



-{Pm\-] 9i{y,t) + QP,g - g{y,t) , (20) 
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Table 1: Zeroth moments of the polarized structure function gi. Also given is the gluon 
component at = 3.0GeV^, AG = 2 /q^ dxg{x). 



M„(MeV) 


Au 


Ad 


As 


AG 


400 


0.64 


-0.14 


-0.01 


0.23 


450 


0.60 


-0.16 


-0.02 


0.21 



The splitting functions Pij are listed in ref |28]. Here it suffices to note that the integrals 



/q^ dzPqq{z) and Jq^ dzPqg{z) vanish at leading order. As a result the zeroth moments 



Aq = 2 I dxg[''\x,t) for q = u,d,s (22) 







do not depend on the momentum scale. For the evolution program being applicable we have 
to assume that g{x, Qq) = which implies that there a no soft gluons at the model scale 
while the effects attributed to hard gluons are approximated by the contact interaction in 
(H). After integrating the differential equations ([T9|)-(pTD from Ql to the singlet and 



non-singlet flavor components are superposed to yield the desired flavor combination. 

6. Numerical Results 

In a first step we determine the model scale Qq. For this purpose we perform the evo- 
lution for the unpolarized structure functions which enter the Gottfried sum rule. As this 
is a non-singlet combination a simplified evolution equation like ( p!U[ ) applies because gluon 
degrees of freedom do not contribute. We vary the lower boundary when integrating the 
evolution equations until maximal agreement with the experimental data available at the 
upper boundary Q'^ = 5GeV^ is obtained. Since details of that calculation will be presented 
elsewhere we content on quoting the result Qg = 0.4GeV^. This is identical to the value 
found for the two fiavor model [jl3l . Naively one would expect that the two and three fiavor 



models would give the same result because the strange degree of freedom cancels in this 
particular combination of structure functions. However, in the two and three fiavor models 
the quark wave-functions are different as can easily be seen from eq (0); in the consider- 
ably simpler two fiavor version the only contributing perturbation is J2a=i{^a/'^)^°', to be 
contrasted with eqs (p!^)- ([l^) . 

All presented results for the polarized structure function at the low momentum scale will 
be taken in the infinite momentum frame (|1^). In table [l| we show the zeroth moments 



(p2D of the polarized structure functions. In ref ||23[ the vacuum contributions to these 
moments were computed. For the sum of valence and vacuum parts those authors give^ 
Am = 0.64, Ad = -0.24 and As = -0.02 using M„ = 423MeV. Apparently the vacuum 
contribution to the zeroth moment is small in accord with our valence quark approximation 
to the polarized structure functions. We also note that due to the deviation from SU(3) 

■^In the present calculation PCAC violating 1/Nc contributions to the non-singlet combinations have 
been ignored. Hence we have to compare our results with the entry "NJL(scalar)" of table X in ref j2^. In 
ref a different regularization scheme to determine the chiral angle and a different expansion scheme of 
the fcrmion determinant were used. This might as well make the (small) differences in the zeroth moments. 
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symmetric baryon wave-functions the quantities shown in table |l| should not be related to 
the F and D parameters determined from semi-leptonic hyperon decaysQ. 

In figure |I| the main result of our calculation is displayed: the strangeness contribution 
g^^"* to the polarized nucleon structure function gi. Apparently the smallness of As does 




Figure 1: The strange quark contribution to the polarized structure function gi for the 
proton in the infinite momentum frame at the model scale = 0.4GeV^ (full line). Two 
cases are displayed, M„ = 400MeV (left panel) and M„ = 450MeV (right panel). Also shown 
is the leading order QCD evolution to = 3GeV^ (dashed line). Furthermore the strange 
quark contribution computed with a pure octet nucleon wave-function (after evolution) is 
displayed (dotted line). 

not necessarily transfer to an overall negligible gi'K A cancellation between positive and 
negative parts is not excluded. This effect is not altered by the DGLAP evolution. We 
furthermore compare to the strange quark contribution obtained with a pure octet nucleon 
wave-function, i.e. with the symmetry breaking effects omitted when diagonalizing (pISl). 
Apparently the incorporation of symmetry breaking effects in the nucleon wave-function 
via exact diagonalization of the collective Hamiltonian yields significantly less pronounced 
strange quark contributions. This is not unexpected in the collective approach |Q. 

In figure ^ we compare the two and three flavor model predictions for the electromagnetic 
flavor combination for which data from SLAG are available at = 3GeV^ [^. We 
recognize that the inclusion of the strangeness degree of freedom yields only minor changes. 
For Mu = 400MeV the results obtained in the two and three flavor models are almost 
indistinguishable. In particular both versions reasonably reproduce the SLAG data pp|. 
The small change of gi{x) when generalizing to flavor SU(3) is unexpected, after all there 
are at least two significant effects associated with this generalization. First, strange quarks 
appear explicitly, second in SU(3) the Glebsch-Gordan coefficients are different and so are the 



^Strictly speaking F and D parameters are only well-defined in a flavor symmetric formulation. 
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g:(x) 



n|=3, Q-=3.0GeV- 
n;=3, Q-=0.4GeV- 
n,=2, Q-=3.0GeV- 
nr=2, Q-=0.4GeV- 




g,(x) 



ii,=3, Q-=3.0GeV- 

ii(=3, Q-=().4GeV- 

ii,=2, Q^=3.0GeV- 

ii,=2, Q^=().4GeV- 





Figure 2: Comparison of the two and three flavor calculation of the polarized structure 
function gi for the proton in the infinite momentum frame. Also shown is the leading order 
QCD evolution to = 3GeV^ Left panel: M„ = 400GeV, right panel: M„ = 450GeV. 
Data are taken from ref ^ . 



nucleon matrix elements of the collective operators^. Apparently these two effects partially 
cancel each other. 

In figure ^ we show our prediction for the gluon component of the polarized structure 
function at the experimental scale = 3GeV^. We remind the reader of the assumption 
that at the model scale {Qq = 0.4GeV^) this component is taken to be zero, i.e. this 
component is solely due to radiation and absorption of soft gluons. We remark that the 
singularity at x — is weaker than 1/x since xg{x) — > in this limit. 

7. Conclusions 

In this letter we have presented some first results for strange quark contributions to 
nucleon structure functions. In particular we have concentrated on the twist-2 piece of 
the polarized structure function gi. On the whole we find the surprising result that the 
contribution of the strange degrees of freedom do not significantly alter the results of the 
pure two flavor model. In both the two and three flavor versions of the NJL chiral soliton 
model we obtain reasonable agreement with experimental data for the polarized nucleon 
structure function gi when the scale dependence is accounted for according to the DGLAP 
scheme. 

To calculate various flavor components of gi at the low momentum scale we have employed 
the valence quark approximation to the NJL chiral soliton model in flavor SU(3). Although 
the valence quark contribution is known to strongly dominate the static axial properties for 
the self-consistent chiral soliton, in the unpolarized case some conceptual requirements like 
positivity are slightly violated by this approximation |jl3|. Hence the full calculation would 



^The inclusion of symmetry breaking effects within the Yabu-Ando approach drives these matrix elements 
towards their SU(2) values. 
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Figure 3: The gluon component of the polarized structure function evolved from the model 
scale Ql to = 3GeV^. The results are displayed for two different values of the up 
constituent quark mass M„. 

require to include the effects attributed to the polarized vacuum as well. As discussed, we 
do not expect them to have significant impact on the result. In addition an extensive next- 
to-leading order DGLAP evolution could be performed for the twist -2 structure functions 
under consideration. As the main effect we would expect a variation of the model scale Q^, 
which anyhow is a free parameter. 
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